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ABSTRACT 


The Approximate Self-Consistent Field Molecular Orbital Theory 
(ASCFT) with complete neglect of differential overlap (CNDO) has been 
applied to the calculation of electronic structure for several selected 
molecules in this text. The theory leads to calculated equilibrium 
geometry configurations and the barrier to internal rotation of ethane, 
inversion barrier of ammonia and the barrier for cis to trans conversion 
for difluorodiazine. These results are in reasonable agreement with 
experimental values in most cases. The stretching force constants are 
higher than experimental values by a factor of approximately two. The 
orbital energy levels for ammonia and for oxygen difluoride as a function 
of bond angles were in reasonable agreement with the full LCAO-SCF 
Calculations. A variation of parameters to find a new set of parameters. 
was attempted with unsuccessful results. 

The Extended Hückel Theory (EHT) can be applied to aliphatic 
hydrocarbons but its application to non-hydrocarbon molecules is quite 


limited. 
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I. INTRODUCTION 


There are many levels of sophistication for the calculation of 
electronic structure using molecular orbital methods. The simplest 
approach is the Hückel Molecular Orbital Theory (HMO). This is one of 
the semiempirical methods which is used extensively in the quantum 
mechanical study of the mobile electrons of conjugated molecules. The 
HMO was extended by many people. One of the theories examined in this 
paper is the "Extended Hückel Theory (EHT)" which was developed by 
Roald Hoffmann (1). A higher level of sophistication is the self- 
consistent field (SCF) molecular orbital calculation. One of those 
methods which has recently been developed is based on a theory intro- 
duced by J. Aw Pople: “Approximate Self-Consistent Field Molecular 
Orbital Theory (ASCFT)" (2), (3), (4). This ASCFT Calculation is 
predominantly used throughout this text. 

We can say molecular orbital methods are "curve fitting". They 
are really based on semiempirical methods. The purposes of the calculation 
described herein are the prediction of potential surface by Pople's 
theory to examine and to try to improve the agreement of a calculated 
potential surface with an experimental potential surface. 

The background describing each of these theories will be discussed 


in the next section. 


II. Httckel Molecular Orbital Theory (HMO) 


In this theory, the mobile electrons can be treated as occupying 
a set of delocalized molecular orbitals (by the Pauli Exclusion Principle, 
not more than two electrons in each molecular orbital). These molecular 
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orbitals are eigenfunctions of a one electron Hamiltonian, The 
orbitals are approximated by a linear combination of atomic orbitals 
(ICAO) centered on various atoms. The Slater type atomic orbitals 
(5) are generally used as approximations to the actual atomic orbitals 
for each atom. When this theory is applied to a wide range of molecules, 
it has the advantage of great simplicity but serious discrepancies result 
because of the neglect of interelectronic repulsions for the case of 
many electrons. 

For the unsaturated molecules: TheJC - molecular orbitals are 
formed by a LCAO of the 2% atomic orbitals which, it is assumed, share 


the same nodal plane. In N-X electron systems, the molecular orbitals 


N 
Wj = 2 Ce Q, | (1.1) 
| (1- птен JM) 


where U represents the jth molecular orbital, Q, represents the atomic 


are given by 


orbital in the de and Gu represents the coefficients of the jth 
molecular orbital and the zth atom (atomic orbital), It is noted that 
the subscript i, j, k, 1 etc. stand for the molecular orbitals and M 22 J 
G, À etc. stand for the atoms or atomic orbitals. 

From a given electron Hamiltonian, Equation (1.1) must satisfy 


the equation 


AU; - ЕД; ES 


Where Ej is the eigenvalue of the jth state, U; is a column matrix, 


and H is the Hamiltonian operator, This can be put in matrix form: 


10 








Y EOS V, 
i O t2 `. O i | 
А | ы; Zu | (1.2) 
U, О--- -о “En Un; 
* 
We multiply equation (1.2) by U: and integrate over all space, dt Б 


then we obtain 


(wa Y AT = (ыш AT (1.3) 


Substituting equation (1.1) into (1.3) 
x 
X GIG Way GOE фт ш 


Definin 
de - 1%, al, dr (1.5) 
and > Ses o, dT (ДЕС) 


We may write КТ (A), 


» e әб» Haz = Е ΕΖ, Ci Суди (1.7) 


Taking οἵ ibi of equation (1.7) with respect to Cer to minimize 


the energy Fi 4 G0 0) , We obtain 


1 Gallo? Gas Προ 8) а 
з = Sap (1.9) 
Haz à Ha (1.10) 


Therefore we obtain the equation 


2, Cg Hu» - Еау Су Sun ΠΠ, 
B L Ca, (Hos » Er = =O 
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51-12) 


This equation (1.12) is the well-known secular equation and gives rise 


to the secular determinants 
| Han m Ek Jal E O (1.13) 


Equation (1.13) is the non-trivial solution of (1.12). 

Now we discuss the Hermitian Property of the Matrices ( H |мл) 
апа Sa A matrix H is said to be Hermitian if it is equal to the 
conjugate of its transpose, i.e., 

H =H 
Where H=(Hu) md H*= (Hx) "x 
VP 2,77 N) 


By definition, a real Hermitian matrix is symmetric, but an imaginary 
symmetric matrix is not necessarily an Hermitian matrix. It is noted 
that the eigenvalues of an Hermitian matrix are all real. 

Next we discuss some transformation properties. In general, the 
atomic orbital ф,, does not form a complete orthonormal set. This 
basis, by » can be transformed into an ο basis by using a 


unitary transformation matrix such that 


new ld 
φ, =<, y ф, (1.15) 


new 
where % is the new orthonormal basis formed by transforming the old 
basis dd in N dimensional space and 2A is the transformation matrix 


element. The old H can also be transformed into a new basis such that 


W 
Πω Š Z ΤΑ! ф 4/2, AT (1.16) 
It may be recalled that the old basis is me. (& H Фф, 47 5 
Therefore equation (1.16) can be rewritten 
Ho = 2, Тае Ре, (1.16)! 
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or 


new * ,,old 7 
На, Е б, Tac Не» bz (1.17) 


GA 


It is noted that if the transformations were made: 


new | 
Фф, = 2, D Tw (41:8 ) 


then the equation (1.17) becomes 
new x. old 
Haz = E қ σα Lu (1.19) 
T) A 
Equations (1.17) and (1.19) can be written by matrix notation such that 


| Ar H "e / (1.37): 
new k ^, , old 
- = | H T (1.19)! 


| / 
From equation (1.19)', if T is a unitary natrix( T -=T >) then we have 


> al (1220) 


This equation (1.20) is frequently used and usually T is real and H91d 


and 


is a symmetric matrix. 


Referring to the secular equation (1.12), we define Age=Ex » then 


2 ( Han Cut m ev" =) (1.21) 


This can be expressed in matrix notation in the following manner 


HC -SC) Га 2 (1461)! 


where H=( Hau ) д! Ark) 
C = (Cuk) 
S « (S5) A 2 N) 


Therefore, we obtain Am 
A -0«&)z2C (S H)C = 
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It is generally true that the matrix (5-1ῃ) may not be symmetric. In 
equation (1.22), the eigenvalue cannot be obtained by methods usually 
employed to diagonalize Hermitian matrices. One way to solve this 
problem is to transform the matrix S to an orthonormal basis, then 
diagonalize the H matrix. The procedures are described as follows: 
Using the transformation matrix Q such that 1р-0ф, we transform the | 


matrix S to the identity matrix I 


S ë Q SEN Te | 
ог Οσο (Q Ώ y i 


Also we diagonalize the matrix S by a unitary transformation matrix T 


such that 


^ - E ا‎ (1.25) 
к сом (FA 
Кта&ту тат) 
Equating eguation (1.21) and (1.26), we obtain x 


ο κι (1.27) 


Then we transform the old basis H to the new basis by the transformation 


|! 


(1.26) 


matrix Q such that 
ER HIE 
-! 0 E 28 
= (TTX) (THAT) -- 


Finally we can diagonalize H'*" by the matrix U. Then we obtain the | 
eigenvalues 


E= UH U’ 
MOTA O NOA ek 


(1.29) 


where U is the matrix which diagonalizes the HNeW and T is the unitary 


matrix which diagonalizes S. Also the correct eigenfunctions are 


W -UTA*TO (1.30) 


The general operations for the application of HMO theory to 
Specific problems are 
(a) Setting up the Secular determinant 
(b) Expanding the Secular determinant 
(c) Finding the eigenvalues 


(d) Determining the M.0.'s by finding the eigenvectors. 


III. Extended Hückel Theory (EHT) 

Hückel calculations originally were applied exclusively for the 
Planar conjugated and aromatic systems of JL electrons with the neglect 
of overlap for non-adjacent atoms. In the extended Hückel theory (EHT), 
Roald Hoffmann (1) expanded the basis set to consist of 2s and 2P 
orbitals for Carbon atoms and ls orbitals for Hydrogen, with the 
inclusion of all overlap integrals. The EHT can calculate molecular 
orbitals for the aliphatic and aromatic compounds which gives good 
results for most hydrocarbons because of the resonable set of parameters 
chosen for the calculations, The parameters are nicely adjusted to the 
self-consistent value for the hydrocarbons but when the EHT is applied 
to non-hydrocarbons, it is not as successful, and even predicts a lack 
of minima in many bond stretching motions. 

III-A. Method of Calculation 


For the calculation of the molecule CnHm, we use a basis set 
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consisting of m-Hydrogen Slater orbitals, n-2s and 3n-2P orbitals. 

The Slater exponents were selected to be 1.0 for Hydrogen and 1.625 
for Carbon throughout the calculations. It is noted that 1.20 was 
selected for Hydrogen in the calculation of the ethane molecule. Those 
values are listed in Table 1. The order of the matrix is In + m which 
gives the number of energy levels of the system and the eigenvalue E's 


for off diagonal matrix elements are retained in the Secular Equation. 


TABLE 1 
Slater exponents for valence Shell S and P atomic orbitals. 
Н 1.0. er Sb or απ 


Li 0.65 C 1.625 О ap 


0.975 1.80 2.45 


1.125 1.95 2.60 
1.30 2.125 


x Present calculations we carried out using 1.2 for Hydrogen. 





The diagonal elements Huu are chosen to be valence state ionization 
potentials; the particular values used are essentially those of Skimner 
and Pritchard (6), For example, the 552 valence state ionization 


potentials for Carbon and ls Hydrogen atoms have the following values: 


Fu (Cap) = — | pa: е. 


Huu (Cas) 
Huu (Hs) = -13.6 еи 


The choice of the off-diagonal term Fz is approximated as 


ALA EV. (2.1) 
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H OS ( Huu + Haz) (Sua) 


The approximation (2.2) was first discussed by Mulliken (7) and used 


(ο ο) 


in the M.O». Calculations by Wolfsberg and Helmholtz (8). 
The value of K in equation (2.2) was chosen to be 1.75. These 


calculations give gairly good results for hydrocarbons. 


IV. Approximate Self-Consistent Field Molecular Orbital Theory (ASCFT) 
Approximate methods for obtaining self-consistent field molecular 
orbitals have recently been developed by J. A. Pople (2), (3), (h). 
Slater type orbitals (5) are used as basis functions. Overlap integrals 
between two different atoms were calculated using the formula given by 
Mulliken, Rieke, Orloff and Orloff (9). It is assumed that the Slater 
type atomic orbitals are real atomic orbitals, and that the M.0.'s for 
the ground state are obtained by Linear Combinations of Slater type 
atomic orbitals. If the Coefficients in the LCAO Molecular orbitals 
are chosen to minimize the total energy, one obtains LCAO - SCF orbitals. 
The best LCAO approximation to SCF functions were first given by Roothaan 
(10). It is noted that the EHT treated all the electrons independently. 
EHI wave functions are used as a starting point for the self-consistent 
field interations in this А5СЕТ, 
The SCF methods were developed using two approximations based on 
the neglect of the differential overlap between all valence atomic 
orbitals. One of these involves complete neglect of differential overlap 
(CNDO) in all basis sets connected by the transformations as mentioned 
in Section I, that is: the off-diagonal elements of overlap matrix 


in the Hamiltonian matrix are to be neglected. The other involves the 
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neglect of the products of qa bu of any two atomic orbitals 
Du and ф, in all electron-repulsion integrals. This complete 
neglect of differential overlap of all basis sets will result in non- 
orthonormal basis sets. 

The CNDO method in ASCFT should be regarded as an approximation 
to a full LCAO - SCF calculation which uses a minimal basis containing 
only inner-shell and valence-shell atomic orbitals. The Slater | 
exponents used in these calculations are listed in Table 1 (for Hydrogen, 
1.2 is used), The values for 3° and ionization potentials in 2s and 
2p atomic orbitals from Lithium to Fluorine were set semiempirically. 
Those values used in these calculations are tabulated in Table 2 and 


Table 3 respectively. 


TABLE 2 


o 
Bonding Parameters /ÓA (їй е) 


Atom Н Li Be B C N О Е 
7 9 m Lf 21 25 EE ο») 


TABLE 3 
"Average" Ionization Potentials Ly, (in. eds) 
H Li Be B C N 
13.06 -- — = -- 
== 9.32 MOS 19... 25.58 32.30 ЫЕ 
5.96 5580-1087 13.19 TOTO 
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The CNDO method of ASCFT can be used for closed ground shell 
structures and doublets, but a triplet structure is not available with 
the computer program used (for example O», By and Co). The computer 
program is presented in Appendix 1. 

IV-1. Self-Consistent Field (SCF) Equation 

The initial discussion covers Roothann's LCAO - SCF method (10) 
as a basis for calculating the MO's of a molecule, paying particular 
attention to transformation properties. The Is Hydrogen and valence 
shell electrons for Li to F are assigned to ICAO molecular orbitals, 
Minimization of total energy by a variation principle with SCF orbitals 
leads to a series of simultaneous non-Linear equations, the so called 
Hartree Fock equations (11), which are applied to LCAO molecular 
orbitals with a closed shell configuration. These equations were shown 
by Roothaan to reduce to the form (10) 


> » €" 5 £ up Coi C 


(4.1) 


where (using atomic units) 


Paz - Bw + Guy (4.2) 
Huv -f 437-2 Waj) h dr (1.3) 
Guo = 2; Р. -z (jic lao) (hel) 


(ард) fed bodas | ше 
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Ju i: J n 2 47 (1.6) 


Do = om Ci (ж (4.7) 


Where Ei is the orbital energy for the ith Μ.Ο. Huz is the matrix 
element of the one electron Hamiltonian including the kinetic energy 
and the Potential energy in the electrostatic field of the core, this 
being written as a sum of potentials Vy (r) for various atoms A in the 
molecule. Gua is the matrix element of the potential due to valence 
electrons and depends on the M.0.'s via the population matrix Parr ; 
The total energy of the molecule is obtained by adding the repulsion 
energy between isolated cores and the separated valence electrons. The 


repulsion term can be approximated by a point charge model 


Erster š electronic 2 2, 2423, (4.8) 


where Ы is the core charge of atom A. The electronic energy is 


given by 


Ba LH +: if -Ка) 


ο] (11.9) 


Buc LA Cui Gil Пир + Cul Gu 


4,2,4 «4,2, ὃ 
px ә Нар *z2 b аи 
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dect Е t2, E Hu» "i Gua | 


11 


72, PualHuo + Bu | (4.20) 


Au 


IV-2. Theory with Complete Neglect of Differential Overlap (CNDO) 
First consider the simplest version of the theory which involves 
the complete neglect of differential overlap between atomic orbitals 


on the same atom. The approximations are given below. 


Approximation A 
The atomic orbitals Pu are treated as if they form an ortho- 
normal set such that Suz = Sur and then the E, forms an 


orthonormal matrix. The condition for orthonormality becomes : 


2, GG 5 924 (шелігі) 


With this approximation, the diagonal matrix elements E 1Α1 (electron 


density) correspond to the electron populations of the atomic orbitals 


ф, and 
2, В, = 2N (4.12) 


where ZN is the number of valence electrons. 


Approximation B. 


Equations (l.l) are neglected because of the orthonormality of 
atomic orbitals. This means (|ж ) is zero unless Y=D and 


А= 6 , The non-zero values will be written Van such that 
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ы = (әл и) 


Approximation С 


The electron-interaction integrals (72 are assumed to depend 


(00803) 


only on the atoms to which the orbitals Du and b, belong and not 

on the actual type of the orbitals. This means that there remains only 
a set of atomic electron-interaction integrals hy measuring an 
average repulsion between an electron in valence atomic orbital on atom 
A and another in a valence orbital on atom B. The set of atomic orbitals 


on one particular atom A is replaced by an alternative set tm, 


E 2 Ory Ф, TRIS 


A general electron-replusion integral in this new basis set is given by 


(tm tn 5.0.) = Z, Oy Oy ( Q Q фф) (11.15) 


Now we apply approximation C, then equation (L. ҚА becomes 








(ъф, da) = {дв отт (4,16) 


It is noted that equation (4.16) is invariant under the orthogonal 
transformation matrix Ó which gives the new basis +. 


Using approximation A, and B, the matrix elements uz become 


Hay = Hau + Ga 
= H+ Be ul) - КАСЫ?) 


(hom 
ΠΝ; AL 


p = Haz - + Data 


(LATÓN 
Equation (1.17) and (h.18) can be rewritten by approximation C 
е а 47 
L. - Huw 2 ا ا‎ Hs Ye (4.19) 


But = Huz - + Pao AB (4.20) 


where ДШ belongs to atom A and Pag is the total valence electron 


density on atom B, 


Να E 2 Ra 


Also the core matrix Mu can be separated as follows: 


Hau (al- 27 Y а) 2 (ul Val) 
Ua “Ἢ @ Ар (1.22) 


(1.21) 


(1 


Н (1.23) 
s= Uw- F ТІКЕ) 

Ж E ВА) | Ci?) 

Where ІШ Y is the diagonal element of JA with respect to the 


one-electron Hamiltonian containing only 
the core of its own atom ( سلا‎ is an atomic quantity measuring the 
energy of the atomic orbital) and Όλο is the one-electron matrix 


element using the local core Hamiltonian ( J| and ZZ are on the 


2з 


same atoms). It is noted that the cores consist of the nucleus and 


inner-shell electron. 


Approximation D 
The interaction terms of the distribution ф, dy on atom A 
and the nucleus of atom B, ( ul Ve |» ) are a to zero if NAD. 
The integrals VIVAP ) are taken to be ho same for all valence 


atomic orbitals on Atom A and there results 


Vas = (ы Vs |) (h.2h) 


It should be noted that the attraction term matrix element Vie 
between nucleus B and the adjacent valence electron /M in the field 
of the nucleus A need not be symmetric. Therefore the equations (h.22) 
and (4.23) can be rewritten 


Fu = سلا‎ - 22,06 V on lm A) 


(4.25) 


На - О πα) ) both Қы» (1.26) 


Same Atóm 


iG as noved that Hu» 4 O if Pu and tu are on different 


atoms which will be discussed in the next approximation. 


Approximation E 
Tr b. and Dy are on different atoms, it is only necessary to 
consider the possible lowering of energy levels by being in the electro- 
static field of two atoms simultaneously which is referred to as a 


"Resonance Integral" and is denoted by Suv $ Huw can be 


2h 





approximated as 


а = Luz = Bae Sun (1.27) 


О 
where „дле is a parameter depending only on the nature of the atoms 
A and B, Empirical values will be used for these parameters. 
The yp matrix elements now reduce to the form from equation 


(1.19) to (h.27) 


kau = Usu (Par ~ 2 Pay) Van FALLE - 


е | 
а = Зав Su». — z- lu SAB (4 
(аза) - 
where ф, belongs to atom A and ἂν belongs to atom B. 


Using the same approximations, the total energy is given by one 


or two atom types from the equation (1,8) 


P total -2 Ea 5 2,2, Ж 
A<B 





Where Т 
кер ыи ee FP) 
1.30) 
de = 22, (2а –2-Раїйв) 
Coe +P Pasa - Pa Ma - Poo Vs) ανν 
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IV-3. Choice of Parameters of CNDO1 and CNDO2 

In the complete neglect of differential overlap approximation, 
the CNDO1 and CNDO2 methods are different levels of approximation 
of the equations. 


The LCAOSCF Hatree-Fock matrix elements are 


ο nz fem Ve) aan 


buy = Sw — ¿Pus ln (ux z) (h.29) 


O 
where the bonding parameter 8348 is given by 


Bas = 0.5 (Ba +88) (4.32) 


In the equation (1.32), JA is chosen empirically but depends only 
on the nature of atom A. These values are listed in Table 2. The 
parameters [25,25 and 2p, 2p can be estimated from ionization 


potentials, using the theoretical value of Ххх (12). 
Is(X,28"2P") = Е(Ху25' 2p") -E(X, 25°2P") 
= — (25,25 — (men-o) Txx (1.33) 
Ip(x,2s"2p") . (x5 2s"2g")- E(X,2s" 2^) 
= — LJop 2p —(Myn-1) Üxx 
аша m Una wit (24-1) 904 (4.35) 


) 
and the atomic electron affinities Au S are used by the formula 


—Au F up + ZA AA (ТЕЗ 
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(1.31) 


or 








Also the equation (1.28) can be rewritten 


» = LJ. + (5ια- ΣΡ) бад 
Ш y (Re- Zp) ав 


BA) 
= DB, (Za дв = VaB) i 
BEA) Q 
where the last term of equation (l.28)! is described as penetration 


integral contributions to huu 
The following table gives the approximate form between CNDOl and 





CNDO2. ED 5 
TABIE li 
Difference between CNDOl and CNDO2 
reference CNDOL CNDO2 
| 
— Un m +(Za-l) бл π. 5) 24-220 x 
(1.37) | 
а Оаа) Ц Ва Bo) a 
ΕΣ (Βεῖρα-- * 7 (Ris Za) (в 
E m Ив) T. (11.38) 


Where 1 the average ionization potentials la (e.v.) used to fix 
Јум are listed in Table 3. 
2 the matrix elements ¿(Is + As) and $(Ip + Ap) are listed 
in Table 5. 


AA TIAS o mte πο 


u asa IA PLE P pates «p AU PII TIP A qr ANTI LO EL ICL APA, I ARTI ST t MP AI a πα 


TABLE 5 





Matrix Elements From Atomic Data (e.v.) 


H πα Be B б N 0 F 
|Қтағда) 7.176 3.16 5,946 9.59 14.051 19.316 25.390 32.272 


Bays) 15058 22563 1,001. CET & 7.275 9.131. 1.688 


em 


If the penetration integral term in equation (.28)' was neglected 
to avoid the excess bonding and high value of force constants, then 
equation (1.28)! can be reduced to the form of (.38). This corresponds 


to substituting the following term for CNDO2 from CNDOL: 


Vas = ZB AB (1.39) 


Using parameters specified in the above way, the LCAOSCF equations are 
solved in the following steps: 
(a) An initial set of M.O. coefficients Cin is obtained from 
a'Hückel type” calculation such that Puu is replaced by 
average ionization potential and Puz by 8 "Su ° 
(Ы) Electrons are assigned in pairs to the M.O.'s with lowest energies. 
(c) Calculated y and Bus are used to form a new Hatree Fock 
matrix element luz) ° 
(d) New coefficients 145. are derived from this | uo and the 
process is repeated from step (b) until self-consistency is 


achieved on all coefficients with a tolerance of 0.0001. 


V. Discussion of Results 

The molecules selected for these calculations are as follows: 
Hydrogen molecule (H5), Lithium molecule (Li2), Beryllium (Be2), 
Nitrogen molecule (N2), Fluorine molecule (F2), water molecule (H50), 
Ammonia molecule (NH3), Ethane molecule (C»Hg), oxygen difluoride 


(Fg0), Difluorodiazine (N5F,), Methane (CH, ), Methylfluoride (CHF) 


and Hydrogen Peroxide (H70,). These calculations will show the minimum 


configurations for bond lengths, bond angles, force constants, the 
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energy level considerations for NH3 and F90, and the Barrier for NH3, 
СоНс, NoFo and Ho05. Finally the calculations will extend the variation 
of parameters to fit the experimental potential surfaces. 
V-1. Comparison of the equilibrium configurations and force constants 
Ihe ASCFT calculations with various inter-nuclear distances and 
bond angles are useful in predicting equilibrium bond lengths and bond 
angles. These results are listed in Table 6. The equilibrium configur- 
ations (bond lengths and bond angles) are in fairly good agreement with 
experimental values but the force constants are too high compared with 
the experimental values by a factor of about 2.0. This can be improved 
by adjusting parameters, such as the Slater exponent and 8° ° 


Sample results are discussed in Section V-9. 


PP SE OE 





| --- -------- 45 „= ° > 


TABLE 6 


Comparison of Bond Lengths, Bond Angles & Force Constants Ar & Re) 
13) ΥΡ kr: 10° Aune Em ) 
Bond Length (aol Bond Angle (deg) (13e ΕΣ. IO" erg rad? 


| Cal. ^ "Obs. бат. OS. | Cal. | (s. 


0.7459 0.7415 | — = 10.143 | 5.6 | 
"uS ЖЕ 


1.1395 1.0976 52.160 m 












"s * Co 


17.260 | 8.2 
104.22 104.45 O.Bhe | Обо ОШ 


м. A ЙН 
ENT 14.68 6.86 | 
Es RER a — 


wa a rit 


5.145 

0.198 15 

5,35 8 
A 

0.63 
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V-2. Water molecules (H50) 

For Ho0, Coy symmetry was assumed with an O-H bond length of 
0.96A°, The first calculation shows an energy minimum at a bond angle 
of 107.209, At this bond angle, the equilibrium O-H bond length of 
1.0289A9 is obtained. With this bond length, the second calculation 
shows the bond angle of 10.22% at which the O-H bond length is 
calculated to be 1.028949. This is in satisfactory agreement with the 
experimental bond angle of 101.150 (13). These calculations demonstrate 
the O-H bond stretching and angle HOH bending motion such that the O-H 
bond lengths tend to have symmetrical stretching motion and the bond 
angle tends to have scissoring motion (symmetrical in-plane deformation). 
À series of these calculations is listed in Table 7. The stretching 
foree constant is calculated to be 17.26 x 10? (dyne/cm) which is 
higher by a factor of 2.1 with respect to the experimental value of 
8.23 x 105 (9упе /ст), The bending force constant is calculated to be 
0.845 x 10-11 (erg/rad?) which is higher by a factor of 1.22 than the 


experimental value of 0.69) x 10-11 (er&/rad2) (11). 


TABLE 7 
Calculation of Roy and HOH for H50 
Calculation I Calculation II Experiment 
Roy (AP) 0.96 (Exp) 1.0280 0.958 
ZHOH 107.28 101.22 101.45 


* The computer outputs (sample) are listed in Appendix 2, 


V-3. Ammonia Molecule (NH3) 
The ammonia molecule with an NH hond length of 1.05A° (13) is 


calculated to be non-planar with an JZ HNH bond angle of 106.825° 
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which is in close agreement with the experimental value of 106.6° (13). 
At this value of the bond angle, the NH bond length was found to be 
1.0657 A9. The subsequent calculation shows a bond angle of 101.869 

and again the bond length of 1.0673 A9. On the other hand the equili- 
brium bond length for the planar configuration gives 1.05h3 A9. It is 
interesting when calculations presented here are compared with results 
recently calculated by R. Moccia and L. Randaccio (17). They reported 
the minimum bond length of 1.878 a.u. (0.9938A°) at a bond angle of 
105%36! and a bond length of 1.837 a.u. (0.9721 A°) at a planar 
configuration. The change in minimum bond length between non-planar 

and planar configuration is 0.0217 A° by Moccia and 0.011) A° by the 
ASCFT. The smaller value of this bond length change in the present 
calculation is due to the higher stretching force constants for ASCFT, 
Also, subsequent calculations can predict the symmetrical bond stretching 
motion and symmetrical in-plane deformation of bond angle for the ground 
state of this molecule. The stretching force constants are calculated 
to be 1h.6h x 10? (dyne /cm). This differs by a factor of 2.1 with 
respect to the experimental value of 6.86 x 10° dyne/cm (U) and the 
bending force constants are calculated to be 0.95 x 10-11 erg/rad.? 


(factor of 1.58 with respect to experimental values of 0.6 x rl 


erg/rad^ (1h)). 
It is also interesting to compare the barriers calculated by the 
several methods and the difference in barriers with respect to adjusting 


the various parameters. These are listed in Tables 8 and 9, respectively. 
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TABLE 8 


Comparison of the resuits of Caiculations for NH3 





алас parrier 
Reference Method of calculations! RNH (qy)| < HNH (de Height(a.u.)Height(a.u. 


Present Calc ЕК 1.9923 |10699! 0.719 0.0251 | 


Higuchi (18) 






Single Determinan | | | | 
Kaplan (19)) polycenter SCF MO |1.916 |1069! | 0.1719 ШО | 


ingle Determinan 

one center SCF MO | 1.928 |110%, | 0.6215 
Single determinant 
one center SCF MO 1,867 [109931 0.626 0.0411 


From Table 9, it is interesting to compare the change in barrier 










Moecia (20) 





Joshi (21) 





with respect to a change in bond angles at the fixed bond length of 
1.0657A° (Calculation II and IV) such that the change in barrier is 
0.000030279 a.u. (0.002510830 - 0.002510551) and the change in bond 
angle is 1.965° (106.825° - 101.86°). From this calculation, it is 
found that the calculated barrier is 0,002510830 a.u. which is different 
by only 1.2% from the first calculated barrier. 

The electronic configuration is found to be 

(18,2 (18) (22,2 bh. 

Configurations calculated by ASCFT, EHT and full LCAOSCFT by H. D. 
Joshi (21) are compared in Fig. 1 and Fig 2. There is complete agreement 
between ASCFT and full LCAOSCFT calculations for the orbital energies 


as a function of height of pyramid. It is noted that the lowest orbital 
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ORBITAL ENERGIES (a.u) 


Fie. 1. Etot. AND Eorb. oF NH3 AS RUNCTION OF PYRAMID HEIGHT. 
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ORBITAL ENERGY (E.V) 


ado 
© 


-200 





Тео. Eist.& Eom.or NH3 vs, PYRAMID HEIGHT 
(BY EXT) 0 


O, | 02 0.3 0.4 0,5 


HEIGHT OF PYRAMID СА”) 


-152.0 


152. 


e.V.) 


TOTAL Ἔν ERGY ( 


کی ووی و ر ر gre‏ = 


energy in LCAOSCFT calculation is due to the (1s) electron (inner- 
shell) of Nitrogen but ASCFT neglects that part. On the other hand, 
the EHT calculation (Fig. 2) does not even demonstrate a double minimum 
potential curve. Therefore the EHT is quite inferior for this kind of 


molecule. 


—sa ο 





Change in Barrier for NH3 w.r.t. minimizing the parameter 


Run СА?) ZZNH | 2НМН Etotal (4.4) | Barrier (asu)! 


Calculation I| 1.02 (fixed) 120 | -13.85761125 | 0.02013903 | 
106.825 | -13.87775028 | 


120 | -13.86130636 | 
(Equilibrium) | 1.0657 106.825 | -13.8891 1215 | 002510551 | 


ο 120 -13. 86,306 36 
= -13.88890154 | 90213 | 
IV 





1.0657 O 128 -13.86362706 
w шый. 101.86 | -13.88971166 | «02608760 
experimentar 1 | 132,206 | 1066 


* 15° (13.88971128 - 13.081308%6) (asus 0.0251,0830 


A A TA A ne pn => =- — sn 
*#— — 9 











V-h. Ethane molecule (C,H,) 
Calculations based on ASCFT and EHT have been carried out for 
Ethane in both eclipsed and staggered configurations with experimental 


values such that Rec E 1.5Ң3А°, Bog = 1.10249 αμα. αθώιΗ = 109937: 


(tetrahedral angle). The first calculation shows the equilibrium angle 
for both eclipsed and staggered forms using EHT and ASCFT. The energy 
then was minimized by varying the bond length Roc» and finally by varying 
bond length Коп. The stretching force constants (Re and сн ) меге 


found and those values are listed in Table 10. 


It is of interest to look at the barrier resulting in the change 
in geometry. These changes in parameters and the corresponding barriers 
are compared in Table 11. These comparisons show that the eclipsed form 
is more open than the staggered form by a few tenths of a degree; also 
the bond length Rog is elongated a few thousandths of an Angstrom and 
the bond length Rog is shortened a few ten thousandths of an Angstrom. 


Some researchers have shown an elongation for Ron. Although these changes 


(Rem) in geometry are small, their effect on the internal-rotation 
barrier is very significant, 
TABLE 10 


Comparison of Ethane Betn. EHT & ASCFT Calculation 




















































. Reference ASOPT ЕНТ Expt. (16) 
clipsed agger clipsed Stagger; ΟΙ 

4. CCH νου 110.86 110.99 110.l | 109937: 
ZHCH 107.60 108. Ol 107.90 108.48 | 1099371 
Rz 1.1660 1.4633 1.8070 1.8000 | 1.543 

RCH 1.1193 1.1195 0.9086 0.9088 | 1.102 
RLCCH Ë (1.74) (1.67) cn «ue |. or63 

e 

лан ЖШС) — ο μ.μ 

Kr ο 15.887 16 . 31,0 197.0 181.0 1.57 

Аса | 54 | 12.125 12.12} 161,05 186.56 | 5.35 
Barrier 2.6586 2.875.5 


(he) 











(23) 





There is some uncertainty in the Table 11 for EHT calculation, 


Every reported calculation scheme has the same tendency to increase 


the barrier when the parameters are changed in sequence, but EHT 
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calculation shows the barrier is decreased. The resultant barrier by 
EHT is 1.6358 X“al/mole. This is not in agreement with the value 
2.6586 Kal mole found by varying only one parameter, ZCCH , holding 
the other two parameters Rec and Rom constant. It is noted that 

R. Hoffmann (1) reported h.O keal mole with a Slater exponent (H) of 
1.0. Present calculations use the value of 1.2. However the results 
are in fairly good agreement with experimental values using both ASCFT 


and EHT (with some uncertainty). 








TABLE 11 








Comparison of Barrier for Ethane 


Parameter (Bel. —Stag) || Barrier (Eclip, —Stag) 
2CCH (deg) | Kee (4) | Ren (AD) ||zccH |] Rece | Ren — | 
L.S. Bartell 0.20 0.0029 0.0001 A ΤΕΝ 
(1960) 
J.T.Vanderslice**%4 0,17 | 0.0024 -0.0001 0.387 | 0.416 | 0.43 
& E.A.Mason(1960) 
27 

МШ” |8 | оо ое | ола | олов | олот 

| (28) с 
x 0.40 0.0059|-0.0002 || 0.779 | 1.108 | 1.298 | 

oye x 

ne 0.55 0.0080| -0.0011 || -1.532 | 1.675 | 2.9kk 
Oe Je Sovers & 30 





Reference 












ien 0.42 0.0027 | -0.0002 1.8637 | 2.2818 | 2.3288 
(ЕНТ) a 0.0070 | -0.0002 2.6506 | 1.5746 | 1.6356 


Expt. 3.03% 0.3(KC@2 /more) (31) 2.875+0.125 (Keal/mole) | 
(23) 
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V-5. Methane (CH),) and Methylfluoride (CHF) 

Calculations between methane and methylfluoride were carried out 
using tetrahedral symmetry for methane and C3v symmetry for methyl- 
fluoride with Вс = 1.11354? for CH), and RSG = 1.109 and сене 5 
Ғог CHF. The calculated values for the equilibrium configurations are 
listed in Table 12. The force constants are calculated for each molecule 
and these are listed in Table 12. The stretching force constants are 
different by a factor of 2.3 for CH), ande? .5 for CH3F with respect to 
their observed values (15). 

For a methane molecule, the calculated bond lengths and bond angles 
are the same as the observed values. 

For methylfluoride molecules, the bond angle ZHCH is more open 
than the angle of < HCF which can be explained as follows: there exists 
rehybridization for the CH4F molecule because there must be more 
P-character in C-F bond due to the electron withdrawing fluorine; on 
the other hand, there must be more S-character in C-H bond which gives 
more SP? hybrid character for the C-H sigma (S) bonds. It is obvious 
that the SP^ hybrid orbital is more planar than SP3, Therefore 2 HCH 
is more open. It is here that the best agreement between theoretical 


and calculated results was found. 
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TABLE 12 


Comparison of f , © and kr ko between CH), € CH3F 








Bond Length (A°) Bond Angie (deg) | 









A MAA A AAA Сеоне a T Ἵκισκαν-ταν 


Cal. Obs. 
Be de cl E а m] 


Ка 1.1180 | 2.109 | _ 1.7115) 


ess [ss] sas 
Οὐ ὅκα |. | 120.382 | 109.997 | 
жа | | 108.5 | 108.98 | (1.77) |0.792(15) 


V-6. Difluorodiazine (NoF>) 








Calculations have been carried out for difluorodiazine in both 
trans and Cis configurations with Кым = J eo Ки = 1.1149 and 
ZNNA = 115° (13). In both configurations, the new equilibrium 
distances and bond angles are determined and are listed in Table 1h. 
When comparisons are made of the equilibrium geometry in Table ly, the 
bond angle 2NNF is more open for Cis than trans configurations. 
Significantly, the bond length Ryy is elongated a few thousanths more 
for trans than Cis (same for RR NF D 

It is of interest to compare the energy level difference between 


Cis and trans configurations. The experimental results giving the 





stability between both configurations are listed in Table 13. In 

1963, Armstrong and Marantz published that the Cis form is more stable 

by 3.0 $ 0,3 Kcal/mole (36) (refer also (37)). The present calculation 
x predicts the Cis form to be more stable by 3.249 *““l/mole than the 


trans (Coh) form. Fortunately, these results are in most satisfactory 


er 


agreement with experimental value. In Table 1) are compared the diff- 
erences in energies between the Cis and trans configurations for a 
variation in molecular parameters, First the angle “ NN was varied 

to minimize the energy. With this bond angle, the minimum bond length 
Ryn was calculated and finally the minimum bond length Ryr was found. Іп 
each case, the difference in energies was compared with the one obtained 
initially. It was found that there is no effect due to varying the 

angle ZNNK and bond length RNN but there is a very significant 
difference when the bond length RNF was minimized. This result shows 


the configuration is affected significantly due to the Ryp stretching 


motion. 






TABLE 13 


Results of experiments for NF, (stability) 

















Author Reference Results 
iie Dif, Bauer (32) J.Am.Chem.Soc. ,69, 3104 trans is more stable 


C.B. Colburn,et.al. "non n" " — ,81,205 than Cis configuration 


R.H. Sanborn ^^ JeChem.Phys., 33,1855 |Cis and N-NXg& is 
"n 


more chemically active 
R, Ettinger,et.al. 













3 
а Р " , 3,2187 lisomer 
G.T. Armstrong 35 .Chem.Phys., 38,169 [Cis is more stable by 
€ Se Marantz 3.020,3 kcal/mole 
— [Cis is more stable 
by 3.249 Xcal note 


The next interesting calculation is the barrier due to the rotational 
motion which was reported by the Binėnboyn (37) . It was reported by 
experiment to be 32.0 keat йыла and calculated to be 8 kcal/mole (36) 
between trans and linear configuration. However, the present calcu- 
lations were carried out in various ways and those results are listed 


in Table 15. The calculated barrier between linear and trans form in 


ho 
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Table 15 is too high with respect to the reported barrier but the 


calculation between trans and C2 form shows a satisfactory agreement 


with the calculated (38) and experimental value (39). 


It must be 


noted that when an attempt was made tocalculate the rotational barrier, 


no dihedral angle exists; that is, when the dihedral angle is 90 degrees, 


the bond angle АМЬ turns out linear, and also the dihedral angle does 


not exist and the motion is predicted to be an up and down motion which 


is shown in Calculation IV and V in Table 15. 


are illustrated in Fig. 3. 


The calculated barriers 


TABLE 15 


Comparison of Barrier for NoE, 


Parameters i | 
No. Ref. Btot. (a.u.) 


п 
| cet) [amor 
Фе (—-—-3|180.0 | 1.209 


! Е(ІТ) - EIN) E(TI) - τπτ) E(I) - EIN) E(I) - E( 
Barrier = 694156 = 201.356 = 72.7057 = 207,651 ^ nct 
S E(II) -- E(V) = 32 kcal/mole (37) 


Cal. (IT) -- E(V) = 8h kcal/mole (36) 


Fig. 3 Barrier for NoFoas function of XNF 
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V-7. Oxygen difluoride (F20) 

For the Fo0 molecule with Cov symmetry (RGF = 1.41A° (13), the 
lowest energy configuration is calculated to be bent with ап ДЕЕ =97.16° 
(at fixed Rop-1.50A9) as compared with the experimental angle of 103.89 
(38). The bond length RS (at angle of 97.169) is calculated to be 
1.179249. Also the bending and stretching force constants are calcu- 
Latedutoube deco ο ο квк aoa UR еге / ad 
and 0.9 x 10? dyne / m (observed RGF =5.57 x 10? dyMe /om (UO). 

It is of interest to compare the orbital energies as a function 
of internuclear angle with that publisted by Buenker and Peyerimhoff 
(39). The calculated results are illustrated in Fig. h and these were 
compared with calculations by Buenker and Peyerimhoff (using full 
LCAOSCF theory) also in Fig. h, The comparison in Fig. h shows a 
satisfactory agreement except for 58}, ΠῚ and 355. It is noted that, 
in full LCAOSCF calculation, the lowest doubly degenerate orbitals 
(due to ls electrons of two Fluorines) and the next singly degenerate 
orbital (due to the ls electron of Oxygen) were neglected. Therefore 
the comparison is valid starting with the third level in the Fig given 


by Buenker. 
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L1G. 4. ORBITAL ENERGY OF 20 AS FUNCTION OF ZFOF. 
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V-8. Hydrogen Peroxide (H>0») 

For the hydrogen peroxide molecule, the calculations were performed 
to find the barrier and the equilibrium configurations with Rgg70.97A?, 
ZH0021009, R2571.h9A9 and dihedral angle of 93.69 (13). Unfortunately, 
the present calculation does not show the existence of a dihedral 
angle. The result of the calculation shows that there exists a Cis 
and trans configuration for H505 just as for N5F5. The structural 


parameters are illustrated in Fig. 5. 


Fig. 5 Structural Parameters (Hz O2) 
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V-9. Homonuclear diatomic molecules 
For the diatomic molecules, the equilibrium configurations of H2, 
Liz, N2 and F2 were calculated. An attempt was made to find a new set 
of optimized parameters by changing the Slater exponents and beta 
zeroes an The equilibrium bond lengths and force constants are 
listed in Table 16. It is noted that B2, Ce and 02 molecules are excluded 


in the calculations because they are ground state triplets, and the Be2 


molecule ig non-bonded. The necessary partial derivatives were found and 
i 
listed in Table 16, 


TABLE 16 






Equilibrium Configurations and Some Partial derivatives 
for homo nuclear diatomic molecules 


ΛΟ Rr(0 do) average | 
[ Obs. | cal. | Obs. [she | ( 


нә 10.716 |0.7ı [10.14 | 5.6 Dm ο 
(0.89) en [^ 
(2.788 (o. P". 


For those molecules, an attempt was made to find a new set of 












= 






optimized parameters fitting experimental equilibrium bond lengths 

and force constants. The Slater exponent and 8 were considered in 
adjusting the parameters to optimize agreement for calculated potential 
surfaces. The procedures used to find a new set of Slater exponents 
and beta zero (,9 ) were as follows: 


First: Change the Slater exponent by # 0.1 and beta zero by £1.0, 


then find the calculated equilibrium bond length and force constant. 


16 


Second: 


Ihird: 


Fourth: 


Find the partial derivatives | 


e ) GF, 85.9, м τ. 


where "S", stands for Slater exponent. 
Find the new Slater exponent and by following equations: 


R= Gas + (6 2,4 (5.1) 
dkr: Ө dS + (95) 546° (5.2) 


then approximate the dR, ds and as follows: 
dR = AR = Reu. — Kal. 


AS = AS = Snew — Sail. 

o Ó 
O s 44 = new – Вел. (5.3) 
Solve equation (5.1) апа (5.2) to find the new Slater exponent 


and yom . For example: 
AS 
58° 
Snew = 0,923 


(1) Ho molecule: 
Bnew = 9.457 (5.4) 





e -5.Ң -3 باه‎ - 7 
-11.83 19.18 1.055 0.157 


therefore 


(2) Lin case: 


-0.49 294} 0.022 -0.16) AS 
-0.59 ( 820ү 0.52 м -0.410 I 46" 
therefore 


Snew = 0.119 
Bnew = 8.59 (5.5) 


(8) N5 case: 


ONB5 «Ομ. -000 0.87 45 
-29.06 L 61.06 2.31) -35.5 (46° 


therefore 
Snew = 2507 
Bnew = 2965 (5.6) 
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From these results, (5.h)and (5.5) are consistent but (5.6) 
is completely inconsistent because the force constant is too 
high. 

Fifth: Finally, repeat the calculation using the above new set of 
input parameters (Snew and ΚΘ nex). The result for the Ho 


and Lio molecules are listed in parentheses in Table 16. 


VI. CONCLUSION 

The Approximate Self-Consistent Field M.O. theory was applied 
to the calculation of electronic structure for several selected mole- 
cules and these results were compared with results by full LCAOSCFT 
calculations, if they were available. The most important features of 
these calculations are to predict equilibrium configurations (bond 
length and bond angle) and the barrier by the variation of internuclear 
distances and angles using the ASCFT, The EHT calculations are limited 
tc hydrocarbons which give good results (refer tocalculation of Ethane 
molecule). In general, the barrier predicted by ASCFT calculation is 
in satisfactory agreement with the experimental value. 

An attempt was made to find a new set of input parameters (Slater 
exponent and «8° ) to improve the prediction of potential surfaces 
for certain homonuclear diatomic molecules but the attempt failed 
except for the Ho and Lio molecules (see Table 16). 

The features of Potential Surfaces of molecules can be predicted 
with the approximate self-consistent. field theory. Qualitatively, the 
agreement between calculations and experimental observations is quite 


good. It is conceivable that, with the adjustment of certain parameters 
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such as resonance integrals (6 9 Slater exponents and penetration 
integrals, quantitative agreement with experimental results can be 


obtained. 
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